Graphene as a reversible spin manipulator of molecular magnets 
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One of the primary objectives in molecular nano-spintronics is to manipulate the spin states of 
organic molecules with a d-electron center, by suitable external means. In this letter, we demonstrate 
by first principles density functional calculations, as well as second order perturbation thoery, that 
a strain induced change of the spin state, from S—l — >■ S=2, takes place for an iron porphyrin (FeP) 
molecule deposited at a divacancy site in a graphene lattice. The process is reversible in a sense that 
the application of tensile or compressive strains in the graphene lattice can stabilize FeP in different 
spin states, each with a unique saturation moment and easy axis orientation. The effect is brought 
about by a change in Fe-N bond length in FeP, which influences the molecular level diagram as well 
as the interaction between the C atoms of the graphene layer and the molecular orbitals of FeP. 

PACS numbers: 75.50.Xx, 81.05.ue, 75.30.Wx 



Fe-porphyrin (FeP) is a potential candidate among 
metalloporphyrin molecules for biomolecular devices to 
be used for spin dependent electronic transport using 
the characteristics of spin switching. A significant ef- 
fort has been made to switch the properties of organic 
molecules by the underlying magnetic substrates. By a 
combined study of element-specific magnetic measure- 
ments and first principles electronic structure calcula- 
tions, it was shown that it is possible to tune the magnetic 
coupling between Fe at the center of FeP with a ferro- 
magnetic substrate like Co or Ni Also, switching 
properties of molecules can be driven by some external 
means like temperature, pressure, light and external lig- 
ands. It should be noted that the change in the spin 
state of the FeP molecule depends on the ligated or non- 
ligated situation Q. A ligated FeP molecule consisting 
of 5 Fe-3d electrons, can have only odd moments {S— 
1/2, 3/2, 5/2) whereas a non-ligated one, with 6 Fe-3d 
electrons, has even spin states (S= 0, 1, 2). Recently, it 
was shown that the spin of metalloporphyrins exchange 
coupled with ferromagnetic films, can be controlled by 
nitric oxide molecule acting as an axial ligand In a 
theoretical work , an organometallic framework based 
on Cr-porphyrin arrays was proposed as a spin filter in 
spintronic devices. 

Graphene has made enormous sensation in the scien- 
tific community due to its exotic multidisciplinary fea- 
tures 0, which have the potential to be exploited in fu- 
ture electronics, gas sensing, batteries, understanding of 
high energy experiments etc. During the synthesis of 
graphene or at a post-synthesis step involving function- 
alization with e.g. an acid treatment Q, a number of 
defects may arise. In sp^ bonded carbon, the divacancy 
defect can easily be introduced in a graphene lattice [9| 
and these are the defects we propose to utilize here. It has 
been shown that the electronic and transport properties 
of graphene are affected significantly due to the presence 



of divacancy defects 13, U | • The chemisorption energies 
at these defect sites are quite large and hence, divacancy 
defects behave as excellent trapping centers for adatoms 
and molecules. 
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Utilizing the novel properties of FeP and graphene, in 
this letter, we propose an alternative method to manip- 
ulate the spin states of a FeP molecule placed at a diva- 
cancy center in a graphene lattice. Our aim of the current 
study is to investigate the chemical and magnetic interac- 
tions of FeP with a defected graphene lattice and to tune 
the magnetic properties of the molecule by modifying its 
electronic and magnetic states by strain engineering of 
graphene. We have also studied the modification of elec- 
tronic states in graphene due to the presence of FeP. To 
this end, we have performed state-of-the-art first prin- 
ciples calculations that show that the strain engineering 
of graphene lattice modifies the spin state of FeP in a 
reversible way. 

We have performed first principles density functional 
calculation using the VASP code IlSj. Plane wave pro- 
jector augmented wave basis was used in the local spin 
density approximation (LSDA) for the exchange correla- 
tion potential. The plane wave cut off energy used was 
400 eV. A 3x3x1 Monkhorst Pack k-points set was used 
for the integration in Brillouin zone. Atoms were relaxed 
until the Hcllman Feynman forces were minimized up to 
0.01 eV/A. For the FeP-graphene composite system, we 
have used a 8x8 supercell of graphene with 21 A of vac- 
uum length perpendicular to the graphene plane in the 
unit cell. To account for the narrow d-states of Fe in 
the FeP molecule, we have used a LSDA-I-U approach 
where the strong Coulomb interaction is added accord- 
ing to the mean field Hubbard U formalism. The values 
of Coulomb parameter U and exchange parameter J used 
are 4 eV and 1 eV respectively as these values reproduce 
correctly the electronic structure and magnetism of FeP 
in the gas phase fl^ . 
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FIG. 1: (Color online) Total energies (shown with respect to 
the ground state) of FeP in the gas phase, plotted as a func- 
tion of fixed spin moments. The charge densities for IfOMO 
and HOMO-1 energy levels are shown for 2 hb and 4 
cases. Red and blue colors indicate spin-up and spin-down 
electron density, respectively. In the FeP molecule, Fe, N, C 
and H atoms are drawn as red, blue, gray and yellow balls, 
respectively. 



Let us first discuss the spin states of FeP in the gas 
phase. We have performed fixed spin moment calcu- 
lations for the FeP molecule and the geometry of the 
molecule is allowed to relax to the ground state. Our 
calculations reveal that apart from the ground state min- 
imum of the total energy with a magnetic moment of 
2 fj,B, another minimum at 0.55 eV higher energy than 
the ground state, exists with a moment of 4 /is, as seen 
in Fig. [T] The energy barrier between the 2 fiB and 3 
fiB magnetic states is found to be 0.81 eV, indicating 
that one has to overcome at least this amount of en- 
ergy to switch the spin state from 2 /is to 4 fiB- We 
will below refer to these states as the intermediate spin 
state (S=l) and high spin state (8=2), respectively. The 
ground state minimum occurs for a Fe-N bond distance 
of 1.96 A whereas the higher energy minimum has a Fe- 
N bond at 2.03 A . Hence, the spin-state is intimately 
coupled to the Fe-N bond distance. We also show the 
magnetization densities corresponding to Highest Occu- 
pied Molecular Orbital (HOMO) and HOMO-1 energy 
levels for the 2 /i^ and 4 jiB cases. For the S=l spin 
state, the dominant contributions to HOMO and HOMO- 
1 levels arise from d^ and dxy orbitals of Fe, respectively, 
whereas dx2^y2 and dz^ orbitals of Fe are relevant for the 
S=2 spin state. So, a clear change of orbital symmetry 
accompanying the change of the spin state is observed. 

In our calculations for the FeP-graphene composite 
system, FeP is initially placed in such a way that the 
center Fe atom of the FeP molecule is on the top of the di- 



vacancy center and the geometry as well as the spin state 
are optimized to achieve the ground state. The calculated 
binding energy (B.E) is 0.28 eV, which demonstrates that 
the system is chemically stable. In the optimized geom- 
etry, the vertical separation between FeP and graphene 
becomes 3.2 A. The local geometry of the molecule re- 
mains more or less the same as in the gas phase, i.e the Fe 
atom almost remains in the molecular plane with a ver- 
tical deviation of 0.05 A from the molecular plane. The 
Fe-N bond distance is found to be 1.96 A (same as in 
the gas phase) with a magnetic moment of 2.18 hb- The 
small difference in magnetic moment of the free molecule 
with an integer moment of 2 /^s is due to the interaction 
of the molecule with the divacancy site, which makes the 
edge C atoms to be spin polarized. 

In Fig. [2j we show the orbital occupancies of Fe-3d 
states in FeP when trapped at the divacancy defect site, 
in the 0% and 1 % strain situations. For free FeP with 
a total magnetic moment 2 /is, the d^^y orbital is com- 
pletely filled, the dz2 is half filled with one spin-down 
electron and the dyz and dxz orbitals are degenerate in 
energy having two spin-down and one spin-up electron. 
In case of FeP trapped at a divacancy defect site in the 
0% strain situation, the dxy orbital becomes half filled 
with one spin-down electron, the degeneracy between dyz 
and dxz is lifted, dyz becomes completely filled whereas 
dxz becomes partially filled with one spin-down electron. 
The most interesting change is found in the dz-2 orbital. 
The unoccupied part of the spin-up dz2 orbital of the 
Fe atom of free FeP, becomes partially occupied when 
adsorbed. This is clearly observed in the magnetization 
density isosurface plot in Fig. 2. The change in the or- 
bital arrangement is due to the interaction between FeP 
and the divacancy defect of graphene leading to a non 
integer magnetic moment of the combined system. 

We now address the possibility to rearrange the or- 
bital occupation of the Fe-d levels and hence to mod- 
ify the spin state of the system by strain engineering 
of graphene with divacancy defects. We have exerted 
tensile strains by increasing the unit cell vectors of the 
unstrained graphene lattice, considering the cases of 1% 
and 2% strained graphene with FeP placed at the top of 
the defect site and have performed geometry optimiza- 
tions. In case of 1% strain, the B.E of the FeP molecule is 
found to be 0.45 eV. Unlike the case of FeP on unstrained 
graphene, FeP on strained graphene goes through a struc- 
tural deformation due to an enhanced interaction with 
the divacancy center. In the unstrained situation, the 
divacancy defect reconstructs to a 5-8-5 structure and 
a clear 7r-band formation is observed (data not shown). 
But once the size of the vacancy site is increased due to 
the tensile strain, the Pz orbitals of the edge C atoms do 
not take part in the formation of a 7r-band. Instead they 
form more localized levels, with unsaturated bonds. The 
structural rearrangement increases the Fe-N bond length 
from 1.96 A to 2.03 A and this leads to a change in spin 
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FIG. 2: (Color online) Magnetization density isosurfaces for 
FeP on (left) 0% and (right) 1% strained graphene. The iso- 
surfaces have been plotted for an energy window of 0.4 eV 
below the Fermi levels in both cases. The spin-up and spin- 
down densities are drawn in red and blue, respectively. The 
energy levels with the d-orbital character of FeP are shown in 
the extreme left and right for 0% and 1% strained graphene, 
respectively. 



state from S—1 to S=2 of FeP. The energy barrier of 0.81 
eV between these two spin states, observed in the gas 
phase, is overcome due to the interaction between FeP 
and the strained graphene lattice. It is important to note 
that the spin crossover behavior of molecules on surfaces 
is hardly realized in practice due to the strong interac- 
tion between the two systems, [l^. On the contrary, in 
this work, we demonstrate that it is the surface-molecule 
interaction that helps in the spin transition. 

It is interesting to investigate hysteresis effects and 
whether a reversible process is possible via a compres- 
sive strain which may cause a transition from the S=2 
state to the S=l state. If possible, this could lead to an 
ultrathin spintronic device having logic states assigned 
as '0' and '1', depending on the spin state (S=l or S=2), 
which may be influenced by strain effects. To test this 
hypothesis, we released the 1% strain and exposed the 
FeP-graphene entity to a 2% compressive strain. For this 
situation, the geometry optimized structure has the Fe- 
N bond length equal to 1.96 A with an accompanying 
transition from the high spin (S=2) to the intermediate 
spin (S=l) state. Hence, a situation of reversible spin 
switching can be achieved by the application of a rela- 
tively small tensile or compressive strain, as illustrated 
also in Fig. ^ 

To analyze the situation further we note that in the 
1% tensile strained condition, the (1^2 orbital becomes 
completely filled (which is significantly different from the 



unstrained case), the dxy orbital is partially filled (one 
spin-down electron) , the dyz and dxz orbitals become de- 
generate in energy and are partially filled (with two spin- 
down electrons). More importantly, the d^^-yi orbital, 
which is completely empty for the ground state of FeP 
in the gas phase as well as on the 0% strained graphene, 
becomes partially occupied, with this orbital being close 
to the Fermi level, as clearly seen in the magnetization 
density isosurface plot in Fig. 2. The orbital occupations 
are similar to those of gas phase FeP in the high spin 
state. This indicates that the interaction between FeP 
and the divacancy in a strained graphene lattice mainly 
overcomes the barrier between S=l and S=2 spin states 
in realizing a high spin ground state of FeP. The applica- 
tion of a compressive strain (2%) brings back FeP to the 
S=l spin state with the orbital arrangement of FeP sim- 
ilar to a 0% strained case. It should also be mentioned 
that the role of a divacancy in spin switching is crucial. 
We have found from our calculations that without this 
defect, the change in spin state of FeP is not possible 
when a strain is applied to graphene. 

Experimentally, the straining of the graphene lattice 
can be realized by straining the substrate on which 
graphene is deposited and the predicted reversible change 
in the spin state may be achieved, which can be moni- 
tored by the x-ray magnetic circular dichroism (XMCD) 
measurements of electronic structure of Fe 3d states [l^l . 
The discussion of the spin moments of FeP measured 
by XMCD requires special attention as the contribution 
of the spin-dipole term (T^) is known to be quite large 
17| for this class of organic molecules deposited on sub- 
strates. Therefore, we have calculated the 7{Tz) values 
following the method prescribed by Wu and Freeman. [3| 
The values are -1.89 and 1.32 for 0% strained 
(S=l) and 1% strained (S=2) graphene lattice respec- 
tively. The absolute values are quite large and are similar 
in magnitude to what have been observed by Stepanow 
et al. for CuPc. One interesting point to note is 

the change in sign of l{Tz) due to the change in the spin 
state and hence, change in the occupation of d-orbitals of 
different symmetries. This dramatic change will strongly 
affect the effective spin moment [2Sz + 7{Tz)) measured 
by XMCD. 

The consequences due to the change in the spin state 
of FeP on graphene under strain are expected to be man- 
ifested in the uniaxial magnetic anisotropy too. To study 
this, we have considered the effect of spin-orbit (SO) in- 
teraction through the hamiltonian Hso ~ C^.S, where 
^ is the SO coupling constant with L and S being the 
orbital and spin angular momenta respectively. Accord- 
ing to second order perturbation theory, the SO contri- 
bution to the energy is : Ego = -C E« o '^"'^'^'"^'^ 



Here |o) and \u) correspond to occupied and unoccupied 
states, respectively, with \u/o) = \lm,a), where I, m and 
a are azimuthal, magnetic and spin quantum numbers. 
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FIG. 3: (Color online) Schematic diagram showing the energy 
levels and their orbital characters for (left) 0% and (right) 
1% strained cases, considered for the calculation of MAE (see 
text). The filled and empty arrows indicate occupied (HOMO 
and HOMO-1) and unoccupied (LUMO and LUMO+1) d- 
states of Fe, respectively. The big arrows indicate the direc- 
tions of magnetization. As are defined in the text. 



respectively. Such an approach was e.g. used in Ref. |l9|. 
We have considered four relevant states in this calcula- 
tion, HOMO-1, HOMO, LUMO, LUMO-f 1, to estimate 
the SO contribution. All other states would cause a too 
large denominator in the expression of the MAE, and 
hence were neglected. We start by defining the follow- 
ing matrix elements: Ai = (dxy, alli.S \d^2 , a) ; A2 — 
{dj;y,a\L.S\dyz,a) = ((i^2_y2, ^| L.S 1^2, a) ; A4 = 
{dx2-y2,(3\ L.S \dyz,ce) as depicted in Fig.[3]with a and /3 
representing spin-up and spin-down states. For the mag- 
netization along the z-axis, we get = ~ S 
where Ai = 0,A2 = 0, A3 = 0, A4 = § 
E„ - Eo. Hence, E^^(^ 



and 5iE — 



When the magnetization is along the x-axis, the spin 
functions of the local coordinate system can be ex- 



pressed in terms of global spinor functions as: x 



t - 



^--{a -t- /3) and 



^ (^u; ^- ciiiu ^ — ^ (a — /?) . With a similar 
approach as above, we obtain Ai = A2 = A3 = 

A4 = and hence, Ef^ = 0- Therefore MAE = 

^2 

^2 with the 001 axis as 



^so 



the easy axis. In a similar fashion, for 1% strained 

a2 

graphene with S=2, we obtain MAE = 7^ — ^fe 

il I 1 where 

4(i5d,.-Bd^2_^2) 4(i5rf„,-i3<,^,)J 2(E^.^-E^^2_^2)' 

we identify HOMO-1 as \dz2,a), HOMO as \dx2^y2,l3), 
LUMO as \dxy,a) and LUMO-f 1 as \dTr,a). One should 
note that LUMO-l-1 is now a doubly degenerate dj^ state, 
so we need to consider five states to calculate the MAE. 

For MAE calculations, ( for Fe in a 2+ state, has been 
taken as 0.05 eV. [l^ The final results are the following: 
for 0% strain (S^l of FeP), MAE ^ -O.4 meV and the 
001 axis is the easy axis, whereas for 1% strain, S=2 and 
the MAE is 2 meV, with the easy axis being in the plane 



of the molecule. This is a quite interesting result as due to 
the change in the spin state, the easy axis switches from 
out-of-plane to in-plane direction with a 5-fold increase 
in the absolute value of MAE. 

In conclusion, we propose here a new direction of 
combining molecular electronics and graphene science to 
achieve novel switching properties. By first principles 
calculations, we have shown that the S=l magnetic state 
of FeP trapped at a divacancy site of a graphene lattice 
can be switched to a S=2 spin state by applying a tensile 
strain in the defected graphene lattice. One may think of 
other alternatives for the substrate, e.g., polymers, which 
can be stretched too. However, the effect discussed here 
relies on a strong binding between the FeP and the sub- 
strate through the defect site, so that any substrate strain 
translates to the FeP molecule. For graphene with diva- 
cancies this is indeed the shown here, and it re- 

mains to be seen if it also holds for other substrates. The 
doubling of the saturation moment represents a big mag- 
netoelastic effect, which is found to be reversible since 
a compressive strain can bring back the S=l spin state 
from a S=2 spin state. The microscopic mechanism of 
the spin switching is the change in Fe-N bond length of 
FeP and the concomittant change in the electronic struc- 
ture and the level arrangement of molecular orbitals that 
hybridize with defect induced states in graphene. We 
also show that the modification of the spin state is asso- 
ciated with a modification of the out-of-plane easy axis 
to an in-plane one, which is detectable, e.g., in XMCD 
measurements. 

We gratefully acknowledge financial support from 
STINT Institutional Grant for Young Researchers, 
Swedish Research Council (VR) and a KOF grant from 
Uppsala University. O.E. is in addition grateful to the 
KAW foundation and the ERG (project 247062 - ASD) 
for support. We also acknowledge SNIC-UPPMAX, 
SNIC-HPC2N and SNIC-NSC centers for the allocation 
of time in high performance supercomputers. 



[2] 
[3] 

[4] 

[5] 



[8] 



For correspondence: Biplab.Sanyal@physics.uu.se 

H. Wende et al., Nat. Mater. 6, 516 (2007). 

M. Bernien et al., Phys. Rev. Lett. 102, 047202 (2009). 

M-S. Liao and S. Scheiner, J. Chem. Phys. 116, 3635 

(2002). 

C. Wackerlin et al, Nat. Commun. DOI: 
10.1038/ncommsl057, (2010); J. Miguel et ai, J. 
Phys. Chem. Lett. 2, 1455 (2011). 
W. J. Cho et al, J. Am. Chem. Soc. 133, 9364 (2011). 
A.K. Geim and K.S. Novoselov, Nature Mater. 6, 183 

(2007) ; A.H. Castro Neto et al, Rev. Mod. Phys. 81, 
109 (2009); A.K. Geim, Science 324, 1530 (2009). 

V. A. Coleman et al, J. Phys. D: Appl. Phys. 41, 062001 

(2008) . 

M. H. Cass et al., Nat. Nanotech. 3, 676 (2008); A. 
Hashimoto et al, Nature 430, 870 (2004). 



5 



[10] S. H. M. Jafri et al., J. Phys. D: Appl. Phys. 43, 045404 

(2010). 

[11] B. Sanyal et al, Phys. Rev. B 79, 113409(2009). 
[12] O. Cretu et al., Phys. Rev. Lett. 105, 196102 (2010). 
[13] G. Kresse and J. Hafner, Phys. Rev. B 47, R558 (1993); 
G. Kresse and J. Furthmiiller. Phys. Rev. B 54, 11169 

(1996). 

[14] P. M. Panchmatia et al, Chem. Phys. 343, 47 (2008). 
[15] P. Gutlich and H. A. Goodwin, Spin crossover in transi- 
tion metal compounds, vol. I-III, Springer- Verlag. 



[16] J. Stohr and H. C. Sicgmann, Magnetism from funda- 
mentals to nanoscale dynamics, Springer scries in Solid 
State Sciences (2006). 

[17] S. Stepanow et al., Phys. Rev. B 82, 014405 (2010). 

[18] R. Wu and A. J. Freeman, Phys. Rev. Lett. 73, 1994 
(1994). 

[19] D-S. Wang, R. Wu and A. J. Freeman, Phys. Rev. B 47, 
14932 (1993). 



